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Holonomic quantum computation (HQC) may not show its full potential in quantum speedup due 
to the prerequisite of a long coherent runtime imposed by the adiabatic condition. Here we show 
that the conventional HQC can be dramatically accelerated by using external control fields, of which 
the effectiveness is exclusively determined by the integral of the control fields in the time domain. 

Remarkably this control scheme can be realized with net zero energy cost and it is fault-tolerant 
against fluctuation and noise, significantly relaxing the experimental constraints. We demonstrate 
how to realize the scheme via decoherence-free subspaces. In this way we unify quantum robustness 
merits of this fault-tolerant control scheme, the conventional HQC and decoherence-free subspace, 
and propose an expedited holonomic quantum computation protocol. 
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Introduction .— As building blocks for quantum com¬ 
puters, the implementation of quantum gates has re¬ 
ceived considerable research efforts over the recent 
years [1]. It has been reported experimentally that num¬ 
bers of pulse-controlled microscopic systems, such as 
solid-state spins [2] and trapped ions [3] , can be hosts for 
implementation of quantum gates. While enormous the¬ 
oretical strategies for conventional quantum gate imple¬ 
mentation have been proposed, there is a revived interest 
in using geometric phases to perform circuit-based quan¬ 
tum computation, termed as holonomic quantum com¬ 
putation (HQC) [4] , which is enabled by the adiabatic 
quantum theorem. The theorem asserts that at any in¬ 
stant a quantum system remains nearby in its instanta¬ 
neous eigenstate of a slow-varying Hamiltonian, specifi¬ 
cally for a cyclic adiabatic process, a geometric phase (the 
Berry’s phase), is acquired over the course of the cycle [5]. 
The geometric phase is exclusively determined by the tra¬ 
jectory of the system in its parameter space and robust 
against local fluctuation [6, 7]. Consequently, a geomet¬ 
ric strategy for implementation of quantum gates permits 
fault-tolerant and robust quantum information process¬ 
ing. Besides inherent resilience in non-Abelian geometric 
phases [8], HQC has an appealing advantage [9-11] in 
utilizing the state-of-art experimental setups due to its 
close relationship to the circuit model [12-14]. A recent 
experiment has implemented a universal set of geometric 
quantum logic gates with diamond nitrogen-vacancy cen¬ 
ters [15], and evidently it will greatly promote research 
endeavour along this line. 

The heart of HQC is the experimental implementa¬ 
tion of the geometric phase acquired in a cyclic adiabatic 
passage. Despite its advantages, the geometric proto¬ 


col itself is challenged with a dilemma. On one hand, 
any HQC algorithm requires a long characteristic run¬ 
time in order to satisfy the adiabatic condition [16]. On 
the other hand, decoherence or leakage accumulated in 
this long runtime gives rise to errors in the HQC pro¬ 
cessing and may eventually destroy the quantumness of 
the system. To get rid of the dilemma, researchers have 
proposed several different protocols. Over a decade ago, 
Wu, Zanardi and Lidar [17] initiated a scheme by embed¬ 
ding HQC into a decoherence-free subspace (DFS). This 
combined HQC-DFS scheme utilizes the virtues of both 
the fault-tolerance of HQC and the robustness of DFS 
against collective dephasing noise based on the symme¬ 
try structure of the interaction between the system and 
its environment. However, the residual individual noise 
remains and ruins the quantum adiabatic passages dur¬ 
ing the long runtime. Later on the HQC-DFS scheme 
was extended by considering the collective dephasing of 
two neighboring physical qubits [18]. Whereas it is more 
feasible experimentally, this scheme has a more stringent 
requirement for the runtime. Recently a non-adiabatic 
HQC-DFS scheme was suggested where the characteris¬ 
tic timescale is reduced by increasing the characteristic 
energy, at the cost of a harsh restriction for the runtime 
equal to the period of the system [19]. However, the fault 
tolerance from adiabaticity therefore becomes obscure. 

In this Letter, we propose a novel and composite strat¬ 
egy to tackle the long runtime issue in the HQC protocols 
via accelerating the adiabatic passage in DFS. We ex¬ 
plain the mechanism and show specifically that the char¬ 
acteristic timescale of the adiabatic process can be vastly 
reduced by means of external field control [20]. Interest¬ 
ingly, it is found that the particular design or shape of 
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a control function, such as regular, random, chaotic and 
even noisy pulse sequences, is not as decisive as it seems 
to be, but only the integral of the control function in 
the time domain plays the crucial role in speeding up the 
adiabatic passage, which greatly relaxes constraints on 
experimental implementation of these control functions. 
Remarkably, we further discover that our Hamiltonians 
in the adiabatic representation are periodical functionals 
of the integral of the control functions, resulting in a net 
zero-energy-cost control scheme - a new mechanism that 
accelerates adiabatic passages with the same effective¬ 
ness. These lead to a new type of fault-tolerance against 
control fluctuations. 

Decoherence-free subspace for qubit gates .— 
Decoherence-free subspace is based on the symmetry 
structure of the system-environment interaction [23-27] . 
Here we briefly recall the method to realize a universal 
set of quantum gates acting on the DFS as firstly pro¬ 
posed in Ref. [17]. To implement a one-qubit quantum 
gate in DFS, we consider a four physical qubit system 
with the Hamiltonian H = 

where = i((Tf < -f = ^(crf - fffO 

are the XY interactions and Dzialoshinski-Moriya terms, 
is Pauli X(Y) matrix acting on the f-th physical 
qubit and m,l = 1,2, 3,4. This Hamiltonian commutes 
with the operator Z = where erf is a Z 

Pauli matrix acting on f-th physical qubit. By setting 
Jf 2 = Ji 2 cos (p(t), where (p(t) is specifically designed 
for HQC, Jf2 = Ji2 sin (p(t), Jfg = J13 and all other 
= 0, the Hamiltonian becomes 

II(t) = J13-R13 + ’^12 [cosip{t)Ri2 - smip{t)R\2] ■ ( 1 ) 


system will remain in the corresponding instantaneous 
eigenstate. Consequently, a Berry’s phase is given when 
the system passes along a closed loop in the Hamiltonian 
parameter space, which is path-independent. Without 
loss of generality, one can consider a case where the sys¬ 
tem is initially at the n-th ground state \En). It follows 
that in the adiabatic regime ifn = where 7 n(t) is 

the Berry’s phase given by 7 „(t) = i /p (i?„(s)li?„(s)) ds. 
Here we emphasize that for dark states with eigenenergy 
En{t) = 0 , its dynamical phase vanishes and the remain¬ 
ing overall phase is a geometric phase. 

Equipped with Eq. (2), we are ready to construct our 
expedited-HQC-DES scheme. To build up a one-qubit 
gate in DES, we consider a cyclic Hamiltonian with pe¬ 
riod of T. We first consider a single qubit phase gate. 
The Hamiltonian i7i(t) is formally given by Eq. (2) re¬ 
garding 9{t) = a sin ip{t) = where a is a dimen¬ 

sionless undetermined coefficient. The two dark states in 
the DFS for Hamiltonian Hift) read as jDo(t)) = JO) and 
\Di(t)) = cos0{t) jl) — sm0(t) j3), respectively. 

In the adiabatic regime, under the unitary evolution 
U{T) = T exp[—z /p^(isi7(s)] where T is time-ordering 
operator, the dark states \Do) and \Di) become 

e*7o(T) li:'i(T)), (3) 

respectively, where 'yj{T) is the Berry’s phase for \Dj), 
j = 0,1. In this manner we achieve a one-qubit phase 
gate by e*T''>(^)lDo(T))(Do( 0 )l + e*T'i(^)lDi(T))(Di( 0 )l. 
Note that \Dj(T)) = jDj(O)). The gate can be expressed 
by a diagonal matrix as diag([e*'''°(^\ The two 

Berry’s phases for dark states are 70 (T) = 0 and 


The bases for DFS have been identified as eigenvec¬ 
tors of Z [17], as spanned by {JO), Jl), 12), j3)}, where 
JO) = JOOOI) and Jl) = JOOIO) constitute the two or¬ 
thonormal states for a logical qubit and j2) = jlOOO) and 
j3) = jOlOO) serve as ancilla. This DFS scheme is robust 
against collective dephasing described by Z 0 H, where 
B is an arbitrary Hermitian bath operator. It is straight¬ 
forwardly proven that in the DFS, the Hamiltonian (1) 
can be rewritten as 

Hi(t)=sin0(t)(ll) (2]+ 12) (11) 

+ cos0(t)(e-*^(‘) 13) (2| + |2) (3|), (2) 

where 9{t) = tan“^(Ji 3 /Ji 2 ). 

Holonomic quantum computation in DFS .— Consider 
a quantum system whose dynamics is governed by a time- 
dependent Hamiltonian Hft) with instantaneous eigen¬ 
vectors \En{t)) and eigenvalues Enft). The wave func¬ 
tion \'4>{t)) satisfies the Schrddinger equation and can 
be formally written as |V^(t)) = \^nit)), 

where (fnit) = — fg En(s)ds is the dynamical phase. If 
the Hamiltonian varies adiabatically and there is a non¬ 
vanishing gap between the interested eigenvalues, the 


= J sin^ 6»(s)^^^ds = 7r[l - Jo(2a)], (4) 

where Jo{x) is a zero order Bessel function of the first 
kind, respectively. 

This technique is also applicable in realization of a 
single ax qubit gate. To build this gate, we imple¬ 
ment the Hamiltonian in the same DFS yet spanned by 
{|+) J-) J2), |3)}, where |±) = (|0) ± |l))/y2. It is 
written as 

H,{t)=sin9{t){\-) (2| + |2)(-|) 

+ cos 9it) |3) (2| + |2) (3|) . (5) 

In this case, the new dark states are |Do(t)) = |+) and 
\Di(t)) = cos9(t) |—) — sin9{t)e~'“f‘^*'> |3), respectively. 
The transformations of dark states under time evolution 
are still described by Eq. (3), and the qubit gate reads, 

. 71 / 2 /^ COS 71/2 -fsin7i/2\ 

y—z sin 71/2 cos 71/2 y ’ ' 

which becomes the ax-g&te when 71 (T) = zr. 
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Now we turn to the two-qubit controlled-phase (C- 
Phase) gate in DFS. Since each logical qubit consists of 
four physical qubits, eight physical qubits are involved in 
implementing a two logical-qubit gate. Let us suppose 
that one can implement the Hamiltonian 


iL 3 (i) = sin0(t)(|l,l) (2,1| + |2,1)(1,1|) + 

cos6»(t) |3,1) (2,1| -h |2,1) (3,1|) . (7) 

The four dark states of the Hamiltonian employed in 
implementing C-Phase gate are given by \Do{t)) = 
| 0 , 0 ), \D^it)) = | 0 , 1 ), \D2{t)) = | 1 , 0 ), iDsit)) = 
cos6{t) |1,1) — sm9{t) |3,1), respectively. 

Over a period T, the Hamiltonian (7) drives these 
states into |Do(0)) -)■ \Do(T)), |L>i(0)) -)■ \Di{T)), 
|L> 2 ( 0 )) ^ \D 2 {T)) and |L» 3 ( 0 )) ^ |D 3 (T)), so 

that the two-qubit gate is diag([l, 1,1, where 

73 (T) = 7 i(T) in Eq. (4). Tuning the free parameter 
a, one can get an arbitrary phase gate at will, for ex¬ 
ample, 73 (T) = TT requires Jo(2a) = 0 at the first root 
a = 1.2024. 

Control scheme. — We now come to the case where the 
Hamiltonian H[t) is not in the adiabatic regime. Our 
scheme is to implement a control c{t) upon the strength 
of the Hamiltonian such that [20, 22] 


H{t) 


1 -I- c{t) 


Hit). 


( 8 ) 


We first show that as long as the control is sufficiently fast 
and strong, the system evolution will behave in the same 
way as that in the adiabatic regime, specifically the wave 
function becomes proportional to an instantaneous 

eigenstate of H(t). It is interesting to note that this con¬ 
trol scheme hardly depends on the details of c{t) but its 
integral in the time domain, and is a new type of fault- 
tolerance against control fluctuations. Consequently, the 
evolution of the corresponding dark states are shown to 
be a qualified workstation for HQC and this induced adi- 
abaticity will be utilized to realized the expedited HQC 
in virtue of a fast modulation over Hamiltonian. 

We emphasize that the results given in Eq. (4) are 
invariant under the transformation (8), which is one of 
key points of our proposal. 

To determine the effectiveness of our control scheme, 

I 


we now introduce a quality factor 

/= (l-^) X |(i7i(0)|C/(T)|77i(0))|, 

where Jyi is the difference between the ideal phase (4) 
and the phase acquired during a finite runtime T. Ac¬ 
cordingly, we have 0 < / < 1 where / = 1 if and only if 
the process simultaneously has ideal adiabaticity and re¬ 
tain the Berry phase predicted by (4). Figure 1 shows / 
as a function of evolution time T (blue curve) in the ab¬ 
sence of control (c(t) = 0), and as a function of average 
noise kick’s strength (c(t)) for T (red dashed curve) that 
is not in the adiabatic domain. 



FIG. 1. (Color online) Quality factor / as a function of the 
time T as shown by the blue curve, where T > 60 roughly 
corresponds to the conventional adiabatic condition for our 
model, and as a function of the average control strength (c(t)) 
for T = 1 in the nonadibabtic domain. The control c(t) is 
modelled as a train of pseudo-periodical square pulses with 
a fixed period 2At and duty cycle 50% (see example of c(t) 
in inset). The amplitude of the control pulses is given by 
J(l—p(l/ 2 —r)), where r G [ 0 , 1 ) is a uniform random number, 
J > 0 is a parameter, and p describes the randomness of the 
control (we use p = 0.5 and At = 0.005). The Berry phase is 
numerically calculated by = {IIi(0)|?7(r)|IIi(0)). Here 
we used a — 0.7605 and 71 = 7 r /2 


Mechanism of the adiabatic speedup .—To understand 
the mechanism of our expedited HQC scheme we ex¬ 
pand the wave function in terms of eigenstates |E„(t)) 
of the Hamiltonians. The matrix elements of the Hamil¬ 
tonians in the adiabatic representation reads, Hmn = 
(En\Em)exp{iJl^il + c{t'))E^nit')dt') [20]. For exam¬ 
ple, the Hamiltonian (2) is 


Hiit) 


/o 

0 

0 


0 


0 


0 


—(/jsim 0 

^ ( 0 -| 0 sin 26 l)e*‘^W 
Vo (0-f(^sin20)e-*^(*) 


\ 


(0 -h |9jsin20)e-*^W (0 -h |(^sin20)e*^W 

— ifCOS^O — (^COS^ 06^*'^*-*^ 

— COS^ 06“^**^^*^ —(fCOS^d j 

I- 


(9) 


where C{t) = J* ds[l + c(s)]. It shows clearly that the Hamiltonian is a functional of the integral (7(1) 
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(or the average of c{t) in the time domain) i.e., H[C], 
meaning that controlled dynamics does not depend on 
the details of c(t) but exclusively depends on the in¬ 
tegral C'(t). Such exclusive dependence also holds for 
our Hamiltonians (2),(5) and (7), and is a unique fea¬ 
ture of our chosen Hamiltonians whose energy differences 
Enm = Em — En are time-independent constants. These 
Hamiltonians, as shown in its adiabatic representation, 
are incidentally equivalent to the Leakage Elimination 
Operators [21]. Hence, the control is fault tolerant in 
the sense that the fluctuation or noise of c{t) hardly 
contributes to C{t). More specially, by considering the 
propagator from t = 0 to t = St, where <C 1 and 
(c(t)) 1/St 1, we can write the propagator as, 

U{St) = Texp (-f /q* Hi{t)dt^ « 1 - i /g* Hi{t)dt. The 

existence of the fast oscillating factor e**^^*^ renders all 
the off-diagonal elements of the propagator vanish and 
then leaves a Berry’s phase to the amplitudes of jUi) and 
two bright eigenstates. Noticeably this factor pushes the 
evolution of system into the adiabatic regime by decou¬ 
pling all the four eigenstates. It clearly illustrates the 
advantage of our control scheme: one needs not to care 
about the exact control function because only the inte¬ 
gral C{t) contributes to adiabaticity. 



J=15 


II n , ^ J—20, noise | 

0.0 0.2 0.4 0.6 0.8 1.0 

At 


FIG. 2. (Color online) Quality factor / as a function of 
the kick length At for the zero-energy-cost control when 
T = 10 in the nonadiabatic domain. We use control function 
c{t) = J(l—p(l/2 — r))(—1) , where r is a random num¬ 

ber. The blue dash-dotted and red dashed curves are for the 
noiseless control with p = 0. The green solid curve represents 
noise control where each point is calculated with 10 random 
noise realizations where p = 0.5. Triangles and squares sig¬ 
nify points satisfying Eq. (10). Here we used a = 0.7605 and 
7i = 7r/2. 

Expedited HQC with net zero energy cost. — On closely 
looking into its pattern, we find that the Hamiltonian (2) 
is eventually a functional of the exponent i. e., 

i7[e*‘^(‘)]. Because of the periodicity of e^^'E) ^ our control 
scheme allows for an interesting case when {c'{t)) = 0, 
where c'(t) has alternating positive and negative values 
such that the net energy cost is zero. We first illus¬ 
trate that the above-discussed positive control c{t) (with 


(c(t)) 0) can be exactly equivalent to zero-energy-cost 

control c'{t), when c{t) = — Ti) with the integral 

C{t), and c'{t) = 7r^-(—l)*d(t —r^) with C'{t). It is easy 
to show that H[e''^'E)]^ = H[e''^ due to the periodic¬ 
ity of and (c(t)) = 27r/Ari (At^ = r^+i — Ti) but 

{c'{t)) = 0 for each two consecutive pulses. The random 
intervals Ar^ are much shorter than T in reality, and ide¬ 
ally the net energy cost of the c'{t) control sequence can 
be considered as zero when Ar^ approaches zero. We can 
also analyze the equivalence for the rectangular pulses 
sequences. Based on the first order of Magnus expan¬ 
sion of U{St) we can justify [30] that if the single pulse 
strength J ^ I, the off-diagonal terms in evolution U{St) 
become zero when 

JAt = 27rn, n = 1,2,3.... (10) 

In Fig. 2 we show the numerical simulation of the quality 
factor / for fixed T = 10 with net zero-energy-cost con¬ 
trol as a function of control pulse length At. We mark 
with triangles and squares when Eq. (10) is satisfied. The 
green solid curve in Fig. 2 shows the zero-energy-cost 
noise control which is more robust against the control 
“kick” length At, while noise positive control has promi¬ 
nent oscillatory dependence on At which requires a more 
accurate choice of At (and/or J) according to Eq. (10). 

Conclusion. —To cope with the long runtime issue in 
implementing adiabatic passages, we have introduced an 
expedited-HQC-DFS control scheme to accelerate the 
conventional HQC. In the adiabatic representation, we 
show explicitly that the integral of the external con¬ 
trol in the time domain typically the control pulse se¬ 
quences, rather than the details of the control func¬ 
tions, exclusively determines the efficiency of speeding- 
up the runtime, such that the scheme is robust against 
the stochastic errors in control. More importantly, we 
further find that the Hamiltonian in the adiabatic rep¬ 
resentation is a periodical functional of the integral of 
the control. The periodicity motivates us to design a 
net zero energy cost strategy for speedup which is also 
robust against control imperfections. These novel re¬ 
sults are confirmed by numerical results. This obser¬ 
vation greatly reduces the experimental constraints in 
generating precisely-shaped pulses and allows us to use 
even random pulse sequences. By combining the fea¬ 
tures of this scheme with a scalable DFS, our expedited 
HQC protocol brings together the four-fold advantages 
of all-geometrical HQC, decoherence-free subspace, zero- 
energy-cost control, and our fault tolerant scheme, a typ¬ 
ical scalable, fast and fault-tolerant architecture. We 
therefore expect that this perfect theoretical protocol be¬ 
comes an experimental practice. 
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